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's  has  also  been  propagated  along  the  ramp  in  the  presence  of  a  neutral  gas  cloud  with  pressure  >  0. 1  Torr. 
Computer  simulation  pfithe  experiment  uses  an  axisymmetxip  particle-in-cell  code  that  includes  self- 
and  external-  B1  and  Br  magnetic  fields  and  the  inductive  Ee  .electric  field.  Calculated  contour  maps 
of  the  particle  and  current  densities  of  the  beam  at  selected  times  are  presented.  The  computer  particle 
density  profile,  the  ratio  vfl/vr,  and  the  magnitude  of  the  diamagnetic  signals  are  in  qualitative  agree¬ 
ment  with  the  observations! 
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DYNAMICS  OF  A  ROTATING  PROTON  LAYER 
ALONG  A  MAGNETIC  RAMP 


I.  INTRODUCTION 

Reversed-field  magnetic  configurations  with  closed  field  lines  generated  by  energetic  particles  have 
potential  application  in  controlled  thermonuclear  fusion1'9  and  in  novel  particle  acceleration10'11 
schemes.  With  respect  to  controlled  thermonuclear  fusion,  such  magnetic  configurations  appear  attrac¬ 
tive  because  the  closed  field  lines  may  considerably  enhance  the  confinement  time  of  a  fusion  plasma 
and  the  gyrating  particles  can  serve  as  an  internal  source  of  energy  for  heating  the  confined  plasma. 

Initially,  the  research  effort  was  focused  on  the  formation  of  relativistic  electron  layers.1213  How¬ 
ever,  field  reversal  generated  by  relativistic  electrons  is  not  of  interest  for  thermonuclear  applications 
because  of  synchrotron  radiation  losses.  To  avoid  the  radiation  losses,  Christofifos'4  proposed  the  use 
of  ions  with  GeV  energy  instead  of  electrons.  In  the  Christofilos  scheme,  as  well  as  a  similar  scheme 
proposed  by  McNally,15  the  ion  beams  are  generated  by  conventional  accelerators  and  have  very  low 
current  (~1  mA).  Therefore,  the  use  of  conventional,  low  current  accelerators  leads  to  long  buildup 
times  and  requires  extremely  low  particle  loss  rates. 

However,  high  current  ion  sources  that  have  been  developed  recently  make  possible  the  formation 
of  field  reversed  ion  rings  in  very  short  times.  The  first  experiments  to  produce  proton  layers  by  injec¬ 
tion  into  a  cusp  used  proton  beams  generated  by  reflex  triodes18  and  magnetically  insulated  diodes.17 
These  beams  contained  only  1014  -  1015  protons  and  produced  diamagnetic  self  fields  that  were  <3%  of 
the  applied  field.  With  improvements  in  reflexing-electron  ion  sources18  19  and  magnetically  insulated 
diodes,20  21  beams  containing  1016  -  1017  protons  per  pulse  have  become  available.  These  advances  in 
beam  generation  made  possible  the  transient  field  reversal  by  a  rotating  proton  layer.22 
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Because  the  number  of  protons  available  in  a  pulse  is  comparable  to  the  number  required  for  field 
reversal,  the  formation  and  trapping  of  the  rotating  layer  must  be  accomplished  with  small  particle 
losses.  In  addition,  the  axial  velocity  of  the  proton  pulse  must  be  reduced  almost  to  zero  to  achieve 
trapping  of  the  protons  within  the  potential  well  of  a  magnetic  mirror.  The  required  axial  slowing  can 
result  from  inductive  coupling  of  the  beam  self-magnetic  field  to  a  resistive  medium23  or  by  propagation 
of  the  layer  along  an  increasing  magnetic  field,  i.e.,  a  magnetic  ramp.  This  paper  reports  on  the  propa¬ 
gation  of  a  field-reversed  rotating  proton  layer  along  a  magnetic  ramp  and  compares  the  experimental 
observations  with  computer  simulations.  In  the  experiment,  the  propagation  is  studied  both  in  vacuum 
and  in  the  presence  of  a  neutral  gas  cloud  injected  into  the  path  of  the  beam. 

Experiments  have  been  performed  in  which  an  annular  beam  containing  2-6  x  1016  protons  of 
energy  — -1-1.5  MeV  passes  through  a  magnetic  cusp  and  propagates  along  a  magnetic  ramp.  Measure¬ 
ments  indicate  that  the  diamagnetic  axial  self  field  of  the  proton  layer  is  as  high  as  1.6  times  the  applied 
magnetic  field.  However  the  amplitude  of  the  diamagnetic  signal  depends  strongly  on  the  strength  of 
the  applied  magnetic  field  and  on  the  self  fields  in  the  Inverse  Reflex  Tetrode  (1RT)  ion  source.24  2S 
Similarly  it  is  found  that  the  loss  of  protons  to  the  wall  and  the  radial  distribution  within  the  rotating 
layer  are  predominantly  influenced  by  the  self  fields  in  the  source  and  in  the  cusp.  Because  of  the 
strong  sell  fields,  the  particle  dynamics  are  considerably  different  from  those  of  single  particles  in  the 
absence  of  significant  self  fields.  For  this  reason,  the  axial  velocity  of  the  rotating  layer  is  higher  than 
expected  for  the  single  particle  case. 

To  better  understand  the  effects  of  the  self  fields,  the  experiment  was  simulated  by  a  particlc-in- 
cell  computer  code.  In  this  code  the  electron  dynamics  were  omitted.  The  effects  of  self  electric  and 
magnetic  fields  in  the  IRT  were  taken  into  account  by  specifying  the  transverse  and  axial  velocity  com¬ 
ponents  as  a  function  of  radius  at  the  cathode  plane,  i.e.,  as  initial  conditions  prior  to  injection  into  the 
cusp  field.  Although  greatly  simplified,  this  model  permitted  the  study  of  the  dynamics  of  the  layer  in 
the  cusp  and  in  the  ramp  fields.  The  simulation  successfully  predicted  most  of  the  experimental  obser¬ 
vations,  for  example,  the  number  and  location  of  the  losses  at  the  wall,  the  radial  distribution,  time-of- 
flight,  and  qualitatively,  the  behavior  of  the  diamagnetic  signal. 
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Experiments  were  also  performed  with  a  neutral  gas  cloud  in  the  ramp.  This  gas  cloud,  air  or  N2, 
was  injected  by  fast  puff  valves.  It  was  found  that  the  beam  ionized  the  gas  and  induced  azimuthal 
plasma  currents  which  magnetically  neutralized  the  layer.  In  contrast  to  experiments23  at  Cornell 
University,  only  a  modest  axial  slowing  of  the  layer  was  observed.  This  is  attributed  to  the  small  ratio 
of  the  azimuthal  velocity  v„  to  the  axial  velocity  vz,  i.e.,  v«/vz  <  1.  The  time  behavior  of  the  diamag¬ 
netic  signal  after  the  passage  of  the  layer  can  be  explained  by  the  excitation  of  magnetosonic  waves.  In 
addition,  several  features  of  the  diamagnetic  signals  of  the  Cornell  experiment  can  also  be  explained  by 
this  hypothesis. 

In  the  most  recent  experiments  it  has  been  shown  that  by  increasing  the  radius  of  the  beam  at 
injection  and  suitably  shaping  the  electrodes  of  the  IRT  source,  rotating  layers  not  strongly  perturbed  by 
self  fields  and  having  vw/v.  —  1.7  can  be  produced.  Such  layers  appear  to  be  suitable  for  subsequent 
trapping  in  a  magnetic  mirror  with  a  time  dependent  gate  field. 

In  this  paper.  Section  11  is  a  description  of  the  experimental  apparatus,  the  magnetic 
configurations  used,  and  a  summary  of  the  diagnostic  methods.  The  experimental  observations  of  the 
propagation  of  the  rotating  pulse  are  presented  and  compared  with  the  computer  simulation  results  in 
Section  III.  The  interaction  of  the  layer  with  neutral  gas  is  discussed  in  Section  IV.  A  summary  discus¬ 
sion  of  the  results  is  given  in  Section  V. 

II.  DESCRIPTION  OF  THE  EXPERIMENTAL  APPARATUS 

A  schematic  of  the  experiment  is  shown  in  Fig.  1.  The  IRT  is  powered  by  the  upgraded  Gamble 
II  generator.  This  apparatus  is  the  same  one  used  in  the  previously  reported  experiment  on  the  tran¬ 
sient  field  reversal  by  a  rotating  proton  layer.22  However  in  the  previous  experiment,  the  magnetic  field 
configuration  consisted  of  a  magnetic  cusp  followed  by  a  75  cm  long  uniform  field  and  a  single  magnetic 
mirror.  In  the  present  experiment,  the  magnetic  configuration  is  comprised  of  a  60-cm-transition-width 
magnetic  cusp  followed  oy  an  increasing  magnetic  ramp.  The  pulsed  electromagnets  that  produce  this 
applied  field  are  multiple-tapped  coils  that  can  produce  several  field  shapes.  Some  of  the  magnetic  field 
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configurations  used  are  shown  in  Fig.  2.  In  this  figure  the  axial  magnetic  field  at  a  radius  of  10  cm  is 
plotted  as  a  function  of  axial  position.  The  8  configuration  is  most  often  used.  The  k  configuration  is  a 
shorter  system  that  has  the  same  field  shape  as  the  8  configuration  so  that  a  higher  field  can  be  obtained 
for  a  given  energy  stored  in  the  capacitor  bank.  The  17  configuration  produces  a  higher  mirror  peak 
than  the  8  configuration.  The  X  configuration  gives  a  higher  gradient.  In  the  previously  reported 
experiments,  the  (3  configuration  was  used. 

The  electromagnets  are  energized  by  a  capacitor  bank  with  a  nominal  energy  store  of  —1  MJ. 
The  field  rise  time  is  typically  5-10  ms  and  is  long  enough  to  penetrate  the  3  mm  thick  stainless  steel 
walls  of  the  vacuum  chamber. 

The  inverse  reflex  tetrode  source  consists  of  a  stainless  steel  screen  cathode  K  (75%  open  area) 
and  an  annular  polyethylene  anode  A  of  typically  18.5  cm  i  d.,  37  cm  o.d.,  and  100  /am  thickness 
mounted  on  a  37  cm  o.d.,  cylindrical  anode  stalk.  At  the  end  of  the  stalk  opposite  the  anode  is  a  3  mm 
thick  stainless  steel  plate.  Electrically  connected  to  the  plate  is  a  movable  stainless  steel  disc  G  located 
23  cm  from  A.  Electrons  emitted  from  the  grounded  cathode  K  are  accelerated  toward  the  semitran¬ 
sparent  anode,  pass  through  it  and  form  a  virtual  cathode  VC  between  A  and  G.  Those  electrons 
transmitted  through  the  VC  reach  the  grid.  The  rest  of  the  electrons  reflex  between  VC  and  K  until 
they  are  absorbed.  The  protons  are  extracted  out  of  the  plasma  that  is  formed  from  the  plastic  anode. 
When  the  applied  voltage  increases  or  remains  constant,  protons  directed  toward  VC  are  unable  to 
reach  G  or  they  reach  it  with  zero  velocity.  Some  of  these  backward  directed  protons  are  reflected  and 
after  passing  through  the  hole  at  the  center  of  the  anode  are  accelerated  and  join  the  remainder  of  the 
beam.  Typically  the  operating  pressure  in  the  system  is  below  0.7  mTorr. 

For  the  present  set  of  experiments,  the  IRT  has  been  modified  so  that  conical  anode  and  cathode 
electrodes  can  be  used  (see  Fig.  I).  These  electrodes  provide  a  radial  component  to  the  applied  electric 
field.  This  radial  field  can  reduce  or  increase  the  inward  radial  component  of  velocity  that  the  protons 
acquire  in  crossing  the  A-K  gap.  Thus,  the  pinching  effect  of  the  azimuthal  self-magnetic  field  B„  and 
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radial  self-electric  field  £rself  (that  is  caused  by  excessive  negative  charge  near  the  axis  within  the  gap) 
can  be  substantially  controlled  for  most  of  the  pulse. 

It  is  desirable  that  the  performance  of  the  IRT  not  depend  strongly  on  the  applied  axial  magnetic 
field  Zf0  (in  the  A-K  gap).  This  has  been  achieved  by  suitably  choosing  the  thickness  of  the  anode  plas¬ 
tic  and  the  correct  A-K  gap  for  a  particular  B0  and  applied  potential.  As  described  previously,19  for  a 
proper  choice  of  applied  voltage,  anode  thickness,  and  applied  magnetic  field,  the  IRT  operates  with  a 
nearly  constant  impedance  Z  during  an  appreciable  fraction  of  the  80  nsec  (FWHM)  long  pulse.  Typical 
waveforms  of  the  voltage  across  the  A-K  gap  (equal  to  the  applied  potential  minus  the  LI  inductive 
contribution)  and  the  total  current  are  shown  in  Fig.  3.  The  voltage  initially  rises  to  —0.5  MV  in  10-15 
nsec  and  then  increases  to  —1.5  MV  within  60  nsec.  The  peak  total  current  is  —0.9  MA,  and  the  peak 
power  of  the  device  is  —1.4  TW.  The  peak  proton  current  is  250  kA,  and  typical  proton  pulses  have  a 
time-averaged  proton  current  lp  of  65-110  kA;  the  time  averaged  total  current  70  is  —0.4-0. 5  MA  The 
proton  generation  efficiency  defined  as  the  ratio  7,,/70  is  between  15  percent  and  25  percent.  The  ion 
generation  efficiency  may  be  much  higher  because  ions  other  than  protons  are  not  monitored. 

The  principal  methods  for  diagnosing  the  beam  parameters  are  nuclear  activation  analysis,26 
prompt  y-ray  measurements,27  magnetic  loops,  and  witness  plates.  The  number  of  protons  Np  at  vari¬ 
ous  radial  and  axial  positions  is  inferred  from  measurement  of  the  residual  radioactivity  induced  in  car¬ 
bon  targets  via  the  resonant  nuclear  reaction  l2C  {p,y)uN  (/?+)l3C.  The  number  of  delayed  y-rays 
emitted  from  the  annihilation  of  the  positrons  is  corrected  for  contributions  from  the  reaction  l2C 
(d,n)uN  </3+)uL'  induced  by  the  natural  isotopic  abundance  of  deuterium  in  polyethylene.  No  correc¬ 
tion  is  made  for  blowoflf,  i.e.,  the  loss  of  UA  nuclei  prior  to  detection.  The  radial  profile  of  the  beam  is 
determined  by  individually  counting  small  segments  of  a  carbon  target  activated  by  the  beam.  Simi¬ 
larly,  the  loss  of  protons  to  the  wall  is  determined  from  arrays  of  activation  targets  located  at  the  wall, 
and  loss  to  the  on-axis  probe  assembly  is  determined  from  carbon  cloth  targets  wrapped  around  the 
shaft  supporting  the  probe. 
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Time-of-flight  At  of  the  rotating  layer  is  obtained  from  the  delay  between  the  voltage  pulse  and 
the  emission  of  prompt  y-rays  when  the  protons  strike  a  teflon  target  and  induce  the  nuclear  reaction 
19 F  (p,ay)l60.  The  prompt  y  rays  are  detected  by  an  NPM-54  photomultiplier  and  fast  scintillator 
located  inside  a  lead  housing  having  a  wall  thickness  in  excess  of  5  cm. 

The  shape  of  the  prompt  y-ray  pulse  also  provides  information  about  the  duration  and  time 
dependence  of  the  proton  pulse.  However,  because  of  the  sensitive  dependence  of  the  thick  target 
yield  of  the  nuclear  reaction  on  the  energy  distribution  of  the  protons,27  only  approximate  information 
about  the  time  dependence  of  the  proton  flux  can  be  inferred. 

The  on-axis  diamagnetic  signals  8  B,  are  monitored  with  3  pickup-loops  spaced  50  cm  apart  and 
mounted  in  a  common  probe  housing  made  of  glass.  The  pickup-loops  consist  of  6  turns  and  are  center 
tapped.  The  wires  are  brought  out  through  a  thin-walled  stainless  steel  tube  within  the  1.5  cm  o.d. 
glass  housing.  A  passive  mixer  or  the  mixer  in  the  vertical  amplifier  of  the  oscilloscope  is  used  to  sub- 
stract  the  signals  from  the  two  halves  of  each  pickup  coil  so  that  the  electrostatic  components  of  the 
signals  cancel  but  the  magnetic  components  are  added.  The  signals  are  passively  RC  integrated.  Typi¬ 
cal  sensitivity  is  —30  Gauss  per  volt  in  a  50  ohm  impedance  load  with  an  integration  time  constant  RC 
—2  p.sec.  The  probes  are  calibrated  using  a  test  fixture  which  precisely  positions  the  probe  within  a 
Helmholtz  coil.  In  addition  to  the  on-axis  probes,  several  probes  are  located  at  various  axial  positions  at 
the  chamber  wall.  These  wall  probes  can  be  positioned  to  measure  either  the  8 B.  or  the  azimuthal  self 
magnetic  fielo  S B„. 

Additional  information  about  the  spatial  location  of  the  proton  layer  is  obtained  using  "witness  tar¬ 
gets"  of  brass  or  plastic.  X-rays  from  electron  Bremsstrahlung  in  the  IRT  are  monitored  with  a  scintilla¬ 
tor  and  photodiode. 
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III.  PROPAGATION  OF  THE  ROTATING  PROTON  PULSE 

(a)  Experiment 

A  series  of  experiments  was  performed  to  determine  the  main  features  of  the  propagation  of  the 
rotating  beam  along  the  magnetic  ramp.  The  objectives  of  these  experiments  were  to  measure  the  ratio 
v„/v.,  the  loss  of  protons  to  the  chamber  wall,  the  strength  of  the  diamagnetic  signal  SB.  and  the  radial 
profile  of  the  beam  at  several  axial  positions  along  the  ramp. 

The  amplitude  and  axial  variation  of  8  B.  depend  sensitively  on  the  injection  angle  4>  which  may 
be  defined  as  tan  1  v,/v.,  where  vf  and  v_.  are  the  mean  values  of  the  radial  and  axial  components  of 
velocity  at  the  cathode.  This  is  demonstrated  by  varying  the  electrode  angle  «|RT,  (see  Fig.  1)  which 
changes  <t> .  A  comparison  of  the  8  B.  signals  for  « |RT  =  0 °  (shot  381 )  and  « |RT  =  9°  (shot  475)  can  be 
made  from  the  data  given  in  Table  1.  Tabulated  for  several  representative  shots  are  the  on-axis  8  B. 
signals  (SB.,.  8  B:1  and  8  fl.3)  monitored  by  the  3-probe  assembly.  The  probes  are  spaced  50  cm  apart 
and  the  8  B:l  probe  is  located  45  cm  from  the  cusp  for  shot  381  and  50  cm  from  the  cusp  for  the  other 
shots  listed  in  the  table.  The  comparison  shows  that  the  8  B.  signals  are  an  order  of  magnitude  smaller 
for  the  «IRT  =  9 °  than  «!RT  =  0.  Also  evident  is  the  drop  in  SB.  that  occurs  as  the  beam  propagates 
further  downstream.  The  corresponding  field  reversal  factor  £  =  (8fi;l/£.,,  8  B.2/ B:2,  8  £.,/£_,)  where 
B.„  is  the  applied  field  at  the  position  of  the  it  th  probe,  also  decreases  along  the  ramp.  For  the 
«,RT=  0°case,C  =  (1.25,  0.25,  0.06);  for  the«,RT=  9°case,£  =  (0.09,  0.015,  0.005). 

The  difference  in  the  strength  of  SB.  for  the  two  values  of  0[RT  can  be  explained  by  a  difference 
in  the  injection  angle  and  subsequent  radial  oscillation.  For  the  0]RT  =  0°case,  the  qErsc,(  +  qv.B0 
pinching  force  acting  on  the  ions  makes  <f>  large  and  negative  and  the  protons  move  inward  to  smaller 
radius  and  acquire  a  large  transverse  momentum.  In  contrast,  when  01RT  =  9°,  the  applied  radial  elec¬ 
tric  field  nearly  compensates  for  the  pinching  forces  of  the  self  fields  in  the  gap,  and  thus  </>  and  the 
radial  proton  excursion  are  small.  If  the  mean  radius  of  the  azimuthal  current  distribution  becomes 
smaller  as  the  protons  move  inward  and  the  azimuthal  current  does  not  decrease,  then  the  diamagnetic 
self-field  increases. 
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To  determine  the  radial  displacement  of  the  beam,  a  carbon  activation  target  is  placed  15  cm  from 
the  cathode.  The  target  consists  of  narrow  annular  bands  which  are  counted  individually.  The  number 
of  protons  per  unit  area  and  the  number  per  radial  band  as  functions  of  radius  are  shown  in  Fig.  4. 
These  graphs  are  the  average  of  four  shots  (two  shots)  for  the  0,RT  =  9°  (0°)  case.  The  peak  voltage  is 
1.3  MV,  BQ  =1.16  kG,  and  the  total  number  of  protons  in  the  pulse  N  is  normalized  to  2.3  x  1016.  An 
estimate  lor  <b  is  obtained  from  tan  1  A/ 1 5,  where  A  is  the  distance  in  cm  from  the  inner  edge  of  the 
anode  to  the  half-maximum  point  for  linear  interpolation  between  the  radial  bands.  For  wirt  -  0  ° 
<t>  =  - 15 °;  for "IKI  =  9 °.  0  >  <b  »  -4° 

As  a  consequence  of  the  large  transverse  momentum  associated  with  a  large  lrf>|,  large  amplitude 
radial  oscillations  of  the  beam  occur  downstream  from  the  cusp.  Because  the  guiding  center  of  the 
orbit  is  well  off  axis,  the  particles  come  close  to  the  probe  assembly  (radius  8  mm)  and  close  to  the 
chamber  wall.  It  is  estimated  that  the  radial  bounce  is  approximately  16  cm.  Therefore,  large  changes 
in  the  radial  position  of  the  beam  can  take  place,  and  the  amplitude  of  SB.  can  vary  greatly  with  axial 
position.  Observation  of  the  radial  protile  of  the  beam  substantiates  this  explanation.  The  profiles  are 
measured  by  nuclear  activation  analysis  of  carbon  targets.  The  results  for  typical  shots  are  illustrated  in 
Fig.  5  where  the  number  per  band  as  a  function  of  radius  at  axial  positions  near  the  probes  are  shown 
for  «|RI  =  0°  and  B„  =  1.5  kG  (8-configuration).  The  beam  is  close  to  the  axis  near  the  SB. |  probe. 
At  the  8 B-2  probe,  the  beam  protiie  is  nearly  flat  and  fills  the  cross-section  of  the  chamber.  The  average 
particle  radius  is  larger  than  at  the  axial  position  of  the  8B.,  probe.  Further  downstream  near  the  third 
probe,  the  profile  is  still  nearly  fiat.  However,  significant  attrition  of  particles  by  losses  to  the  wall  has 
occurred.  Therefore,  the  small  amplitude  of  SB.,  can  be  explained  by  the  loss  of  particles  and  the 
radial  position  of  the  particles.  Because  of  the  radial  bouncing  of  the  beam,  the  ratio  f  =  8B.i/8B.| 
depends  on  B,t.  This  is  seen  in  Fig.  6  where  £  is  plotted  as  a  function  of  B„  At  low  field.  £  is  large;  at 
high  field.  £  is  small. 

The  loss  of  particles  to  the  chamber  wall  does  not  appear  to  be  large  in  the  lirst  meier  after  the 
cusp.  Typically,  the  losses  are  substantially  less  than  20%  in  this  region  Further  downstream,  how- 
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ever,  as  the  diameter  of  the  conical  chamber  decreases,  the  losses  can  be  high.  These  losses  depend  on 
the  strength  of  the  magnetic  field  because  the  mean  radial  position  and  the  amplitude  of  the  radial 
bouncing  are  smaller  in  a  high  field.  This  is  illustrated  in  Fig.  7  where  the  losses  for  shot  529 
(fi0  =  0.99  kG)  and  shot  542  (50  =  1.7  kG)  are  shown.  These  losses  are  determined  by  activation 
analysis  ot  10  cm  x  10  cm  carbon  targets  placed  along  the  chamber  wall. 

For  the  lower  field  case,  it  is  evident  from  the  large  losses  between  1.0  m  and  1.5  m  from  the 
cusp  that  much  of  the  beam  is  near  the  chamber  wall.  In  addition,  activation  analysis  of  carbon  cloth 
targets  wrapped  around  the  on-axis  probe  assembly  shows  that  10%  of  the  beam  strikes  the  probe  hous¬ 
ing  between  the  first  and  third  probe  when  the  first  probe  is  located  50  cm  from  the  cusp. 

The  effect  of  the  cusp  transition  width  A^  on  8/3.  and  the  wall  losses  was  investigated.  In  the 
normal  8-configuration,  A-  r~;  60  cm.  This  transition  width  can  be  reduced  to  =25  cm  by  inserting  a 
pair  of  coils  straddling  the  cusp  inside  the  vacuum  chamber.  As  shown  by  the  comparison  of  the  data 
for  shots  477  and  480  (see  Table  II.  there  is  little  effect  on  SB..  Similarly  there  is  no  significant 
difference  in  the  losses  to  the  wall.  This  insensitivity  to  Ac,  is  probably  the  result  of  the  overwhelming 
influence  of  the  self-fields  on  the  motion  in  the  cusp.  For  shot  480,  (,  =  1.67.  Thus,  gross  distortion 
of  the  cusp  field  occurs. 

A  scries  of  experiments  was  performed  to  maximize  the  ratio  v„/v..  Since  the  available  stored 
energy  of  the  capacitor  bank  was  insufficient  to  produce  the  desired  magnetic  field  for  the  high  injection 
energy,  the  potential  applied  to  the  anode  of  the  IRT  (which  could  be  as  high  as  1.5  MV)  was  reduced 
to  0.7  -  0.9  MV.  From  the  time  delay  A-  of  the  half-amplitude  points  of  the  leading  edge  of  the  mag¬ 
netic  probe  signals,  it  has  been  determined  that  the  axial  energy  downstream  from  the  cusp  is  less  than 
50%  of  the  injection  energy.  Since  the  energy  invested  to  build-up  the  self  fields  is  small  in  comparison 
with  the  kinetic  energy,  it  is  reasonable  to  assume  that  the  transverse  energy  is  comparable  to  the  I  *  •  n  ii  i  - 
tudinal  energy.  The  maximum  value  obtained  for  the  ratio  v/v.  is  1.5  and  is  achieved  in  the  8 
configuration  with  =  1.7  kG  for  an  injection  energy  of  0.9  MeV. 
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Tests  with  the  steeper  gradient  k -configuration  of  the  applied  magnetic  field  did  not  result  in  an 
appreciable  change  in  vL / v.  as  determined  from  measuring  A r  for  the  leading  edge.  However,  the 
width  of  the  85.  signals  increases  with  axial  position  because  of  the  difference  in  v.  at  the  head  and  the 
tail  of  the  beam.  Specifically,  in  the  k -configuration  with  B0  =  2.9  kG,  f>B:]  lasts  —65  nsec  and  6 B:1 
has  a  plateau  on  the  trailing  edge  and  a  pulse  width  of  —116  nsec  (FWHM).  The  Ar  between  the 
half-maximum  points  of  the  trailing  edge  indicates  that  v.  —  0.45  cm/nsec  corresponding  to  an  energy 
of  0.1  MeV.  Thus,  it  is  likely  that  a  substantial  fraction  of  the  protons  in  the  pulse  have  v/v.  greater 
than  the  value  obtained  for  the  particles  at  the  leading  edge  of  the  pulse. 

To  further  increase  the  ratio  vt/v.,  the  IRT  was  modified  to  produce  a  larger  radius  beam. 
Because  the  radial  magnetic  field  B,  at  the  cusp  increases  proportionally  to  r,  the  azimuthal  force 
experienced  by  a  charged  particle  passing  through  the  cusp  is  proportional  to  r.  Therefore,  higher  rota¬ 
tional  velocities  can  be  obtained  by  increasing  the  radius  of  the  beam.  The  anode  of  the  new  source  is 
a  48  cm  i  d..  58  cm  o.d.,  150  //m  thick  polyethylene  ring.  In  a  limited  number  of  shots  taken  with 
wirt  =  0°.  the  source  typically  produces  —  I  0-1.5  x  10lh  protons  per  pulse.  From  activation  analysis  of 
carbon  targets,  it  has  been  estimated  that  l«f>l  <  12°.  A  new  magnetic  field  configuration,  A.  is 
employed  to  provide  a  steep  gradient  and  high  magnetic  rirror  peak. 

With  the  large  radius  beam,  a  substantially  reduced  v.  is  observed  downstream  from  the  cusp. 
For  8, |=  1.27  kG  (A -configuration)  and  F„  —  0.98  MV,  the  At  between  the  half-maximum  points  on 
the  leading  edge  of  8  B:l  (at  40  cm  from  the  cusp)  and  8  B.j  is  —70  ns.  This  corresponds  to  v.  —  0.71 
cm/nsec  from  which  an  energy  of  —0.27  MeV  is  inferred.  At  B{)  =  1.6  kG,  and  F0  —  1.08  MV,  the 
At  is  nearly  the  same.  In  this  case  V(/v.  —  1.7. 

The  results  of  this  series  of  experiments  unequivocally  demonstrate  that  the  ratio  vL/v.  can  be 
increased  to  the  desired  level  provided  that  the  appropriate  magnetic  field  is  used. 

Especially  interesting  is  the  shape  of  the  8 B.  signals.  The  85.  j  signal  for  B0  =  1.27  kG  is  narrow, 
(—60  nsec  FWHM)  and  has  a  second  peak  delayed  by  —120  nsec  from  the  first  and  about  one-third  of 
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its  amplitude.  For  B0  =  1.6  kG,  the  second  peak  of  Sfl,,  is  almost  as  large  as  the  first  and  is  delayed  by 
—40  ns.  The  6B:i  signal  has  a  width  of  —120  nsec  (FWHM).  For  both  values  of  B0,  &B:2  has  a  single 
peak,  but  its  amplitude  is  —40%  smaller  for  the  higher  field.  A  possible  explanation  is  that  the  second 
peak  in  the  8fl.|  signal  is  the  result  of  reflected  particles.  If  the  reflection  point  for  the  B0  =  1.27  kG 
case  is  near  the  axial  position  of  85.2,  then  no  second  peak  would  appear  in  8 B:1.  Also  observed  is  that 
the  interval  between  the  peaks  is  about  twice  the  At  measured  between  bB;\  and  8 B,2.  At  the  second 
probe  (90  cm  from  the  cusp)  Bl2  —  6.6  kG.  in  the  higher  field  case,  the  same  field  strength  occurs  at 
67  cm  from  the  cusp.  Therefore,  based  on  the  limited  data  set  available,  it  appears  that  the  shape  of 
the  8BZ  signal  is  consistent  with  reflection  of  a  large  fraction  of  the  particles  in  the  beam. 

(b)  Computer  Simulation 

The  dynamics  of  a  proton  layer  launched  into  the  applied  magnetic  field  configuration  used  in  the 
experiment  are  investigated  by  computer  simulation.  This  study  uses  a  particle-in-cell  code  which  has 
been  described  previously28'30  and  is  similar  to  several  other  codes  which  have  been  used  to  study 
field-reversed  configurations.31 04  The  code  has  cylindrical  axisymmetry  (r,  --geometry),  and  includes  the 
self-  and  external-  B,  and  Br  magnetic  field  components  and  the  inductive  E„  electric  field.  These  fields 
are  derived  from  the  azimuthal  component  of  the  vector  potential  A0  (r.z.i)  which  is  calculated  from 
the  self-consistent  azimuthal  current  density  by  solving  the  equation, 

V2A„  —  A  Jr2  —  -  (j.0  J.  (1) 

The  protons  are  space  charge  neutralized  by  background  or  accompanying  electrons,  and  there  are  no 

axial  or  radial  currents. 

Solution  for  A„  is  obtained  on  a  32  x  128  grid  by  Fourier  analysis  in  the  ^-direction  and  by  Gaus¬ 
sian  elimination  in  r.  At  the  boundaries,  A„  is  set  equal  to  zero  because  the  lime  scales  of  interest  are 
shorter  than  the  time  for  the  magnetic  field  to  penetrate  the  chamber  wall.  The  size  of  a  cell  is  adju¬ 
stable  providing  that  A./A,  <  8  for  numerical  accuracy,  where  A.)rl  is  the  length  of  the  cell  in  the  :lr) 
direction. 
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The  proton  layer  is  represented  by  2000-10000  macro-particles  which  are  injected  into  the 
specified  external  B  with  assigned  canonical  angular  momentum  P„.  The  invariance  of  Pv  is  used  to 
calculate  the  azimuthal  velocity  vH  at  a  given  particle  position.  The  total  particle  energy  (kinetic  and 
field)  is  monitored  and  energy  conservation  to  within  1%  is  found. 

The  initial  conditions  for  the  proton  beam  at  injection  are  chosen  to  model  the  beam  at  the 
cathode  of  the  1RT.  The  protons  are  launched  with  a  specified  energy  and  a  given  injection  angle  (Mr), 
that  is  computed  from 

r  —  r„ 

<t>(r)  =  - 4>m.  (2) 

r2  ~  ra 

where  '$>,„{  =  \6(r:)\  is  the  maximum  value  of|<Mr)jand  <t>(rt])=*  0.  The  protons  that  are  emitted  toward 
the  virtual  cathode  and  which  reflect  before  reaching  the  grid  and  then  pass  through  the  hole  in  the 
anode  can  be  included  in  the  beam  by  selecting  an  anode  inner  radius  r,  that  is  less  than  the  true  inner 
radius  of  the  anode  plastic.  If  these  protons  have  slightly  divergent  velocities,  then  rn  >  r,  can  be 
chosen. 

Both  the  proton  current  and  the  potential  applied  to  the  anode  are  approximated  by  triangular 
pulses  (see  Fig.  8).  The  potential  rises  from  an  initial  value  of  400  kV  to  a  peak  value  V0  and  then  falls 
to  a  terminal  potential  of  400  kV.  Both  waveforms  have  a  pulse  width  of  60  nsec. 

At  various  times  during  the  simulation  run  the  code  computes  the  mean  energy,  "fluid"  velocity 
components  and  contour  maps  of  important  quantities  such  as  particle  density,  current  density,  .4(J°1'11, 
-4“"  and  .  In  addition,  the  field  lines  are  plotted,  and  the  particle  flux  as  a  function  of  radius  is 
obtained  for  a  few  axial  positions.  This  allows  a  comparison  to  be  made  between  the  computer  simula¬ 
tion  and  the  experiment  of  the  location  of  the  particle  dump  on  the  wall,  the  beam  radial  profile,  the 
diamagnetic  signal  on  axis  and  at  the  wall. 

Although  the  code  is  generally  in  good  agreement  with  the  experimental  observations,  it  should 
be  noted  that  several  aspects  of  the  experimental  conditions  have  been  ignored  in  the  simulation. 
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Deluding  electrostatic  fields,  the  inductive  axial  electric  field,  and  the  electron  dynamics.  Also  missing 
from  the  code  are  the  converging  chamber  walls  and  the  removal  from  the  system  of  particles  that 
would  strike  the  wall  in  the  experiment. 

The  radial  profile  of  the  proton  beam  predicted  by  the  simulation  is  shown  in  Fig.  9.  The  1.2 
MeV,  120  kA,  proton  pulse  is  produced  from  an  anode  with  r2  =  18.5  cm  and  r,  =  4.5  cm  =  r0  with 
</>„,  =  -10°  and  -15°.  The  anode  is  located  30  cm  upstream  from  the  cusp,  and  the  pulse  propagates  in 
the  8-configuration  with  fl, >=  1.4  kG.  The  position  of  the  experimental  chamber  wall  is  shown, 
although  the  code  assumes  the  wall  is  at  r  =  33  cm  for  all  axial  positions. 

As  expected,  at  15  cm  before  the  cusp,  the  outer  edge  of  the  proton  beam  is  displaced  inward  by 
/  tan  </>„,.  where  /  is  the  distance  between  the  target  and  the  anode.  After  the  cusp,  however,  it  is  seen 
that  <b  has  a  strong  influence  on  the  radial  position  of  the  protons  and  on  the  loss  rate  to  the  wall.  At 
50  cm  from  the  cusp,  the  profile  of  the  15°  case  is  broader  but  also  has  a  smaller  inner  radius.  The 
profiles  for  both  cases  are  peaked  further  from  the  axis  in  the  downstream  targets.  The  protons  in  the 
15°  case  move  outward  so  far  that  between  1.0  and  1.4  m  from  the  cusp  that  they  hit  the  wall. 

A  comparison  of  a  simulation  for  <6„,  =  -15°  and  several  experimental  shots  is  shown  in  Fig.  10. 
For  both  computer  results  and  experiment,  the  number  of  protons  per  unit  area  is  plotted  as  a  function 
of  radius.  In  the  simulation.  2.25  x  I0lh  protons  are  emitted  from  an  anode  with  r:  =  16.5  cm,  r,  =  2 
cm.  and  r„  =  3  cm.  The  1.2  MeV,  120  kA,  proton  pulse  is  injected  into  the  magnetic  configuration 
with  R, i  =  1.4  kG.  At  15  cm  distance  before  the  cusp,  the  computed  profile  is  in  good  agreement  with 
the  data  of  shot  597  for  which  .V,  =*  2.5  x  I0,(’  protons.  Downstream  from  the  cusp,  the  comparison  is 
made  with  the  available  data  from  shots  with  a  smaller  ,Vp.  To  make  the  comparison  with  the  simula¬ 
tion.  the  data  from  these  shots  (#602,  604,  607)  have  been  normalized  to  Sp  =  2.5  x  10"’.  In  Fig.  10, 
the  normalized  data  is  plotted,  but  the  actual  Sp  is  also  given  for  each  case.  Although  the  coarse  reso¬ 
lution  of  the  experimental  points  makes  the  agreement  in  amplitude  less  precise,  the  width  predicted  by 
the  simulation  at  50  cm  from  the  cusp  is  approximately  the  same  as  that  observed  in  the  experiment. 
Further  downstream  the  agreement  is  good  at  small  radius  but  breaks  down  at  large  radius.  However,  it 
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is  clear  from  the  simulation  that  significant  loss  to  the  wall  occurs  between  100  cm  and  141  cm  from 
the  cusp.  The  integrated  particle  flux  which  strikes  the  wall  of  the  vacuum  chamber  agrees  with  the 
results  of  activation  analysis  of  wall  targets. 

Both  the  shape  and  amplitude  of  the  on-axis  diamagnetic  signals  predicted  by  the  simulation  are  in 
fair  agreement  with  those  observed  experiment  'v  at  various  axial  positions.  For  example,  the  simula¬ 
tion  gives  Sfl.,  —  1.5  kG,  8S;2  —  0.5  kG,  and  8fl.,  ~  0.35  kG.  The  experimentally  observed  values 
are  &B:l  —  2.0  kG,  BB:2  —  0.2  kG,  and  8 B:i  —  0.06  kG.  The  low  value  of  SB:]  predicted  by  the  com¬ 
puter  is  apparently  the  result  of  the  oversimplified  model  of  the  initial  conditions.  At  the  position  of 
the  first  probe,  the  radial  distribution  of  azimuthal  current  and  thus  SB.  are  rapidly  varying  functions  of 
axial  position.  Small  errors  in  the  initial  conditions  can  cause  a  large  discrepancy  in  8 B:l.  In  contrast, 
the  predicted  values  of  8 B:2  and  8 B.}  are  too  high.  This  is  largely  due  to  the  fact  that  the  code  does 
not  remove  particles  from  the  beam  once  they  strike  the  wall  of  the  vacuum  chamber. 

Tnc  radial  bouncing  motion  of  the  beam  is  evident  in  the  contour  plots  of  the  azimuthal  current 
density  J„(r.z )  in  the  r-z  plane.  An  example  is  given  in  Fig.  11  for  a  beam  with  V(l  =  0.85  MV,  lp  —  90 
kA,  and  <bm  =  -15°.  Contours  of  J„  are  shown  in  snapshot  frames.  The  first  frame  is  at  60  nsec  after 
the  beginning  of  emission  of  the  proton  pulse,  and  succeeding  frames  are  spaced  by  60  nsec  intervals. 
The  protons  arc  emitted  from  an  anode  with  r2  =  18.5  cm,  and  r{)  =  4.5  cm.  The  chamber  wall  radius 
is  33  cm  and  the  position  of  the  conical  chamber  wall  used  in  the  experiment  is  shown  as  a  reference. 

In  the  figure,  an  arrow  indicates  the  location(s)  of  maximum  current  density,  and  its  magnitude  is 
given  in  each  frame.  The  current  density  decreases  by  20%  of  the  maximum  value  between  successive 
contours  in  directions  away  from  the  peak  current  density.  Since  five  contours  are  plotted  in  each 
frame,  there  is  an  order  of  magnitude  difference  in  current  density  from  the  highest  to  lowest  contour. 

In  (he  60  nsec  frame,  the  entire  pulse  has  been  emiited.  Most  of  ihe  particles  are  moving  inward. 
The  outer  edge  of  the  beam  has  become  concave,  and  the  density  is  peaked  near  the  head  of  the  beam. 
The  peak  JH  is  near  a  radius  of  9  cm.  and  the  average  particle  radius  is  7.1  cm.  At  120  nsec,  the  head 
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of  the  beam  has  spread  outward  radially,  and  has  very  steep  gradients  near  the  leading  edge.  A  long 
wake  has  developed  at  the  tail.  The  head  extends  along  the  radial  direction  from  4.5  cm  to  18  cm,  and 
the  average  radius  is  12.5  cm.  The  mean  axial  velocity  is  1.1  cm/nsec  and  the  average  vs/v,  ~  0.26. 
In  the  third  frame,  at  180  nsec,  the  head  of  the  beam  fills  the  entire  cross  section  of  the  chamber  (aver¬ 
age  particle  radius  is  11.2  cm)  and  particles  are  intercepted  by  the  wall.  In  addition,  the  wake  of  slower 
particles  in  the  tail  is  lengthening.  In  the  240  nsec  frame,  an  appreciable  fraction  of  the  beam  has 
bumped  into  the  wall. 

A  high  magnetic  field  and  smaller  injection  angle  reduce  the  amplitude  of  the  radial  oscillation  of 
the  beam.  Results  of  a  simulation  using  the  8-configuration  with  B0  =  4  kG  are  shown  in  Fig.  12. 
Plotted  are  J„  contours  in  the  r-z  plane  for  a  beam  with  V0  =  0.85  MV,  lp  =  90  kA,  and  4>m  =  —8.6°. 
Key  features  include  the  formation  of  a  steep  gradient  at  the  head  of  the  pulse  and  a  long  wake  at  the 
tail  In  the  third  and  fourth  frames  (180  nsec  and  240  nsec)  reflected  particles  are  seen,  and  the  axial 
motion  of  the  dense  part  of  the  pulse  is  slow.  The  average  v.  is  ~7.3  mm/nsec  and  vH/vr  —  1.1  at  180 
nsec  In  agreement  with  the  experimental  results,  the  simulation  predicts  small  wall  losses  for  these 
injection  conditions 

'•  iie  influence  of  self  fields  on  the  propagation  velocity  is  most  apparent  for  simulations  of  higher 
current  beams.  As  reported  earlier,28  the  self  field  produces  a  somersaulting,  i.e.,  particles  at  the  head 
of  the  beam  are  decelerated  and  particles  at  the  tail  are  accelerated.  The  result  is  a  spreading  of  the  dis¬ 
tribution  of  particle  energy.  Another  effect  of  the  self  field  is  to  reinforce  the  cusp  at  the  outer  radius 
and  to  weaken  it  at  the  inner  radius  of  the  beam.  Particles  at  the  inner  edge  of  the  beam  thus  do  not 
rotate  as  much  as  a  single  particle  would  in  the  absence  of  the  self  field.  At  the  outside  of  the  beam, 
the  particles  see  a  stronger  cusp  and  gain  a  greater  azimuthal  velocity. 

These  effects  are  demonstrated  in  the  results  of  a  simulation  of  a  350  kA  beam  (V  =  1,3  MV, 
<bm  =  -15°,  fl, |  =  I  kG,  8-configuration).  The  distortion  of  this  cusp  is  apparent  in  Fig.  13,  where  the 
field  lines  of  the  combined  self  and  applied  fields  are  displayed  in  the  r-z  plane.  Reconnection  occurs 
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before  the  wake  of  slower  particles  has  passed  the  midplane  of  the  cusp.  Shown  in  Fig.  14  are  contours 
of  proton  density  at  various  times  after  the  proton  emission  begins.  The  peak  density  region  and 
highest  gradient  are  near  the  head  of  the  beam.  A  long  wake  of  slower  p articles  forms  quickly,  and  a 
large  radial  expansion  of  the  beam  envelope  occurs  between  1.0  m  and  1.9  m  downstream  from  the 
cusp.  The  mean  ratio  vN/v_.  increases  along  the  ramp  to  ~0.3  at  2  m,  and  the  average  v.  slows  from 
1.36  cm/nsec  at  20  cm  after  the  cusp  to  1.17  cm/nsec  at  1.95  m  downstream  from  the  cusp.  However, 
the  head  of  the  pulse  actually  accelerates  along  the  ramp.  Just  past  the  exit  of  the  cusp  v.  =  1.16 
cm/nsec  for  the  head  of  the  beam.  Between  50  cm  and  90  cm  from  the  cusp,  the  leading  edge  averages 
v.  —  1.6  cm/nsec,  and  between  90  cm  and  1.75  m,  the  head  moves  with  v.  —  1.7  cm/nsec.  Thus  the 
motion  of  the  leading  edge  corresponds  to  an  energy  about  10%  greater  than  the  peak  applied  potential. 
In  contrast,  the  average  total  velocity  —1.3  cm/nsec  corresponds  to  an  energy  of  —1.0  MV  that  is  less 
than  the  average  potential  applied  to  the  IRT  anode. 

The  larger  radius  beam  propagated  in  the  A  configuration  has  also  been  simulated.  A  simulation 
was  performed  with  Vp  =  I  MV,  lp  =  50  kA,  1.0  x  |0lf’  protons/pulse.  =  2.3  kG.  r:  =  29  cm,  rt  = 
24  cm,  and  <bm  =  -10°.  The  current  density  contours  are  shown  in  Fig.  15.  Large  inward  radial 
motion  is  observed  in  the  first  100  nsec.  The  average  radius  for  the  two  frames  shown  is  15.3  cm  (100 
nsec)  and  13.9  cm  (200  nsec).  At  100  nsec,  the  average  v.  is  —6.4  mm/nsec  and  v./'v,  —  1.6.  A  sub¬ 
stantial  fraction  of  the  beam  reflects  between  50  cm  and  75  cm  from  the  cusp  so  that  at  200  nsec  the 
average  v.  is  negative  (-0.5  mm/nsec).  Also  the  temporal  shape  of  the  computed  fi B:]  and  hB.:  agree 
with  the  experiment,  although  both  are  somewhat  higher.  The  simulation  predicts  65.,  =  660  G  at  50 
cm  from  the  cusp  and  f>B.:  —  180  G.  For  Bu  =  1.6  kG.  the  value  for  fi#rl  is  about  three  times  the 
observed  value,  and  the  predicted  bB:2  is  about  twice  that  of  the  experiment.  The  discrepancy  between 
experiment  and  simulation  is  probably  the  result  of  differences  in  the  initial  and  boundary  conditions 
and  perhaps  the  neglect  of  electron  dynamics.  In  particular,  the  value  of  f>B.  is  extremely  sensitive  to 
the  radial  position  of  the  beam  which  in  turn  is  sensitive  to  the  injection  angle  and  applied  field. 
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IV.  INTERACTION  WITH  NEUTRAL  GAS 

The  propagation  of  the  rotating  layer  along  the  ramp  in  the  presence  of  a  low  density  neutral  gas 
background  was  investigated.  The  purpose  of  these  experiments  was  to  study  the  formation  of  an 
azimuthal  plasma  return  current  and  to  test  the  coupling  between  the  motion  of  the  beam  and  the 
plasma  currents.  The  experiments  with  neutral  gas  are  performed  by  injecting  a  puff  of  neutral  gas  <N2 
or  air)  into  the  conical  portion  of  the  vacuum  chamber.  Fast  opening  valves  of  the  magnelically- 
driven-hummer  type  arc  used.’6  The  valves  are  mounted  on  the  chamber  1.5  m  from  the  cusp.  After 
the  valves  are  opened,  the  gas  cloud  expands  at  approximately  the  sonic  speed.  When  the  front  of  the 
cloud  reaches  a  desired  axial  position,  the  potential  pulse  is  applied  to  the  1RT  and  the  beam  is 
launched.  The  pressure  is  monitored  by  an  electron  receiving  lube  (type  5702)  operated  as  a  fast  ioni¬ 
zation  gauge.'7  Pressure  front  risetimes  are  typically  less  than  1  msec.  The  main  diagnostic  measure¬ 
ments  of  the  beam  propagation  are  observations  of  the  diamagnetic  signals  and  the  prompt  y-ray  pulse 
from  a  teflon  target  located  200  cm  downstream  from  the  cusp. 

When  the  pressure  of  neutral  gas  is  greater  than  100  mTorr.  azimuthal  plasma  return  currents  are 
excited  that  almost  completely  neutralize  the  azimuthal  proton  current  This  is  illustrated  in  Fig.  lb 
where  on-axis  8  Ji.  signals  with  and  without  gas  are  shown.  In  addition,  a  small  axial  slowing  of  the 
layer  also  occurs  when  it  internets  with  gas.  The  time  of  flight  to  a  Teflon  target  located  200  cm  from 
the  cusp  is  (determined  from  the  half  maximum  point  on  the  leading  edge  of  the  prompt  y-rav  signals) 
—  140  nsec  without  gas  and  175  nsec  with  gas  For  the  shot  without  gas,  the  inferred  average  axial 
velocity  is  1.64  cm/nsec  which  is  nearly  equal  to  the  initial  velocity  of  particles  with  the  peak  injected 
energy  -~1.4  McV  l  or  the  shot  with  gas,  the  average  v  is  —1.3  cm/nsec  corresponding  to  about  0.9 
MeV 

A  large  radial  expansion  of  the  beam  follows  its  azimuthal  current  neutralization.  As  a  result,  a 
substantial  number  of  particles  strike  the  wall.  This  dump  onto  the  walls  begins  within  a  hall  meter  ol 
the  front  of  the  gas  cloud. 
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It  is  apparent  from  the  results  of  Figs.  16  and  17  that  in  the  presence  of  a  background  gas  cloud, 
the  SB.  signals  exhibit  oscillation  and  change  of  sign.  In  Fig.  17,  the  8fl:1  is  large  because  the  neutral 
gas  arrives  at  the  position  of  the  probe  after  the  proton  pulse.  The  6flz)  signal  has  a  long  tail  that  is 
attributed  to  the  clamping  of  the  field  by  the  high  conductivity  plasma.  The  period  of  oscillation  is 
longer  for  8 Bzl  (1.4  ^usec)  than  for  8 B;i  (450  nsec). 


The  observed  periods  are  consistent  with  the  excitation  of  magnetosonic  waves.38  The  theory39 
predicts  that  the  period  of  this  mode  is 


2a 


I  + 


(w  c) 


(L(opl)2 


2  a 
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where  is  the  Alfven  speed,  a  is  the  beam  radius,  L  is  a  characteristic  length  and  ay  is  the  ion 
plasma  frequency.  The  plasma  density  np  can  be  estimated  from  2nn  I',,  or,  where  n0  is  the 
background  density,  is  the  average  proton  flux  of  the  beam,  <r  is  the  ionization  cross  section  —1.2 
x  10  16  cm2,  and  r  is  the  pulse  duration  —50  nsec.  The  factor  of  2  accounts  for  the  ionization  by  the 
co-moving  electrons.  Typically  np  —  2.5  x  10'  Vcm3  for  a  background  pressure  of  —100  mTorr.  For 
shots  521  and  51 1,  the  values  of  the  external  field  at  the  position  of  the  probes  are  B:]  =  2.4  kG,  B,2  = 
3.5  kG,  and  B:i  =  5.9  kG.  The  estimated  values  are  2.8  x  I07,  4.1  x  107,  and  6.9  x  109  cm/sec 
and  a  -  10,  18,  and  17  cm  respectively.  Equation  (3)  predicts  rm  -  0.7,  0.9  and  0,5  /usee.  In  view  of 
the  imprecision  in  the  estimate  of  np.  the  values  of  rm  are  in  approximate  agreement  with  the  experi¬ 
mental  values. 


An  upper  bound  can  be  calculated  for  the  electron-ion  collision  frequency  iy  from  the  theoreti¬ 
cal40  41  stopping  length  /  given  by 


/as  I  On. 
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where  np/  tth  —  10,  v./v„  —  4,  and  v.  =  1.2  x  )09  cm/sec.  From  the  experimental  observation  that  the 
beam  has  lost  about  1/3  of  the  energy  associated  with  its  axial  motion  in  a  distance  of  200  cm.  /  >  400 
cm.  This  length  implies  iy  <  I09sec  '.  Critical  clamping  of  the  wave  requires 
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where  cu)H  is  the  lower  hybrid  frequency,  w ^  is  the  electron  plasma  frequency  and  k  =  rr/a  is  the 
wavenumber.  If  n„  =  2  x  lO'Vcm-1,  and  a  =  10  cm,  Eq.  (5)  gives  the  •—  1.5  x  1010  B  (B  in  kG). 
Because  the  value  of  vtl  determined  from  the  stopping  distance  is  much  less  than  this  value,  the  mode 
is  underdamped,  as  it  is  observed  experimentally. 


Recently  reported  results  of  the  ion  ring  experiment  at  Cornell  University  include  the  observa¬ 
tions  that  ( 1 )  the  azimuthal  plasma  return  current  is  absent  in  H2  or  D:  background  gas  but  it  is  present 
with  He,  Ni,  or  air  background  gas,  and  (2)  the  diamagnetic  signals  have  oscillatory  tails  that  change 
sign.42  These  results  are  also  consistent  with  the  excitation  of  magnetosonic  oscillation.  For  a  nitrogen 
plasma  of  density  —1  x  I014,  B.  —  8  kG,  and  a  =  13  cm  the  predicted  period  of  the  oscillations  is 
—0.43  /xsec.  This  value  is  within  a  factor  of  1.5  of  the  observed  period.  For  a  hydrogen  gas  back¬ 
ground,  the  plasma  density  is  probably  less  than  for  N2  by  about  a  factor  of  3.  Therefore  in  hydrogen 
gas,  for  B.  =  8  kG,  the  period  is  —66  nsec.  Because  this  period  is  less  than  the  pulse  length,  the 
plasma  return  current  is  expected  to  be  small  as  observed  experimentally. 

V.  SUMMARY 

The  experimental  observations  show  that  the  amplitude  of  the  radial  oscillations  downstream  from 
the  cusp  depends  strongly  upon  the  initial  conditions  and,  in  particular,  the  injection  angle  <t>  These 
oscillations,  in  turn,  influence  8  B.  and  the  loss  of  particles  to  the  wall.  Injection  into  the  cusp  at  large 
radius  and  with  a  high  field  and  steep  ramp  result  in  a  ratio  v„/v.  —  1.7  and  in  the  reflection  of  a  sub¬ 
stantial  fraction  of  the  pulse. 


Computer  simulation  of  the  experiment  also  indicates  the  effect  of  <b  on  the  radial  bouncing. 
Good  agreement  is  obtained  between  computed  and  observed  radial  profiles  at  various  axial  positions. 
The  code  predicts  the  axial  position  and  magnitude  of  the  losses  to  the  wall,  and  the  predicted  ampli¬ 
tude  of  HB.  is  in  approximate  agreement  with  the  experimental  measurements.  Simulation  for  high 
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field  and  moderate  <t>  reveals  low  losses  to  the  walls,  reflection  of  some  of  the  particles  as  well  as  the 
large  v„/v_.  observed  for  large  radius  injection. 

Furthermore,  during  the  interaction  of  the  layer  with  neutral  gas  at  moderate  pressure,  the  azimu¬ 
thal  beam  current  is  neutralized  by  the  beam  excited  plasma  current.  The  observed  oscillatory  behavior 
of  the  diamagnetic  signals  can  be  explained  by  assuming  the  excitation  of  the  magnetosonic  mode. 

In  conclusion,  the  main  features  of  the  dynamics  of  a  field-reversed  rotating  proton  layer  along  a 
magnetic  ramp  are  reasonably  well  understood.  It  appears  that  suitable  magnetic  field  configurations 
can  sufficiently  slow  the  axial  motion  of  the  layer  to  permit  the  formation  of  a  stationary,  field-reversed 
proton  ring  that  can  be  used  to  confine  plasma. 
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Table  I  —  Data  for  Selected  Shots 


me.* 

(kG) 

SB., 

(kG) 

Sfl.s 

(kG) 

fttf., 

(kG) 

Shot 

ft 

1.16 

0.15 

0.04 

0.02 

475a 

ft 

1.27 

2.18 

0.69 

0.09 

381 

ft 

1.02 

2.18 

1.01 

0.32 

477 

8  + 

1.05 

2.56 

1.28 

— 

480 

V 

1 .32 

1.66 

0.19 

— 

557 

V 

1.20 

1.92 

0.12 

— 

585 

V 

1.19 

2.24 

0.38 

0.19 

558 

V 

1.14 

1.64 

0.15 

— 

563 

ft 

0.69 

2.18 

1.76 

0.58 

541b 

s* 

0.76 

2.21 

1.26 

0.45 

484b 

s* 

1.05 

2.18 

1.35 

0.75 

482b 

8* 

1.46 

2.27 

0.15 

0.06 

485b 

ft 

1.7 

1.4 

0.23 

- 

542b 

‘m  e  =  magnetic  configuration 

ft  cusp  width  50  cm 

ft 1  cusp  width  ■-  25  cm 

a  H mi  —  d'  lolher  shots  —  O'*) 

h  anode  i  d  =  0  5  cm  (other  shots  =  18  cm) 
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GOLDI  N,  MARSH,  PASOUR.  PLRSHING  AND  KAPLT.AN.AKOS 


POSITION  (cm)  -*■ 


Tig  4  —  Number  nf  protons  per  unu  urea  and  the  number  per  radial  band  striking  a  segmented  carbon  target  15  cm  from 
the  cathode  I  15  cm  before  the  cuspl  with  ■  I  1  b  kC  i  The  two  eases  are  for  a  planar  I K  I  ni|ltl  -  O')  and  fora  conical 
IRT  b'lai  -■  9  I  The  radial  extent  of  the  anode  is  indicated  for  each  case 


GOLDEN,  MARSH,  PASOUR.  PERSHING  AND  KAPETANAKOS 


0.0  0.5  1.0  1.5 

B0  (  kG )  at  the  IRT 


I  ig  h  —  Ratio  o!'  nn-axis  diamagnetic  field  IlKI  cm  front  cusp  lo  Dial  51)  cm  from  cusp  for  tico  field  configurations 
[he  lines  are  least  square  tils,  anil  the  anode  inner  radius  r,  is  gicett  l**r  each  case 


M) 


Number  of  Protons  (10 


70  nsec  t  =  60  nsec  t  =  50  nsec  t  =  0  nsec 


GOLDEN,  MARSH.  PASOUR.  PERSHING  AND  K  A  PET  A  N  A  KOS 


NRL  Ion  Ring  Experiment 

Computer  simulation  of  a  l  =  350  kA,  Ep  =  1.2  MeV  proton  pulse  passing 
through  a  magnetic  cusp. 


7. 

Etg.  13  -  Lines  of  total  magnetic  field  <#„,  +  at  various  times  alter  proton  emission,  as  calculated  he 
simulation  with  V„  -  I  3  MV.  -  350  kA  15  >  I  k(i  In  configuration  I 
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in  units  <>l 


RADI 


cusp 

AXIAL  POSITION  (M) 


lag  I*  -  l*n>tt>n  density  contours  from  simulation  using  a  large  diameter  IR1  (r2-  29  cm.  r,  -  24  cm)  and  a  large,  steeply 
ramped  magnetic  field  I  Hu  2  J  k(i.  \  -conliguration >  Here  V2  =  1  MV.  In  -  50  kA.  and</>„,  =  10  Arrows  locale  maximum 
densitv  \Xl.  and  each  contour  represents  a  change  of  0  2  \^r 
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Shot  493 
No  Gas 


100  nsec/d  iv 
640  G/div 

750  G/div 


2  *  50  cm  after  cusp 
B0  =  1.65  KG 

100  cm 
2.4  KG 


430  G/div 


150  cm 
4.0  KG 


Shot  521 

N2  Gas  (P~  lOOmTorr) 

200  nsec/div 
640  G/div 

750  G/div 

215  G/div 


fig  16  —  Diamagnetic  signals  on  axis  from  a  shoi  without  gas  ami  a  shot  with  gas.  The  top  trace  in  each  set  has  its  polarity 
res cr sell  at  the  scope  lor  display  purposes  In  shot  521.  the  from  of  the  gas  cloud  is  located  -  30  cm  downstream  from  the 
cusp 
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8B  (  Gauss  ) 


GOLDEN.  MARSH,  PASOUR.  PERSHING  AND  RAPE!  ANAKOS 


DIAMAGNETIC  SIGNALS  WHEN 
GAS  FRONT  AT  Z  =  45  CM 


215  + 


150  cm 
B0  =  5.9  kG 


200  nsec/div 


Shot  511 

Ng  (  P~  200  mTorr  ) 

I  ik>  n  —  Diamagnetic  signals  on  axis  from  a  shoi  in  which  the  gas  Iron!  <  200  mlon  air) 

has  just  reached  the  first  probe  <  50  cm  from  cusp) 
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